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Classical Langevin dynamics of a charged particle moving on a 
sphere and diamagnetism: A surprise 
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Abstract. - It is generally known that the orbital diamagnetism of a classical system of charged 
particles in thermal equilibrium is identically zero - the Bohr-van Leeuwen theorem. Physically, 
this null result derives from the exact cancellation of the orbital diamagnetic moment associated 
with the complete cyclotron orbits of the charged particles by the paramagnetic moment sub- 
tended by the incomplete orbits skipping the boundary in the opposite sense. Motivated by this 
crucial, but subtle role of the boundary, we have simulated here the case of a finite but unbounded 
system, namely that of a charged particle moving on the surface of a sphere in the presence of 
an externally applied uniform magnetic field. Following a real space-time approach based on the 
classical Langevin equation, we have computed the orbital magnetic moment which now indeed 
turns out to be non-zero, and has the diamagnetic sign. To the best of our knowledge, this is 
the first report of the possibility of finite classical diamagnetism in principle, and it is due to the 
avoided cancellation. 



In this work we re-visit the problem of the absence of 
classical diamagnetism of a system of charged particles in 
thermal equilibrium. This vanishing of the classical dia- 
magnetism in equilibrium is generally referred to as the 
Bohr- van Leeuwen theorem [1-4] . The fact that classically 
the orbital diamagnetic moment vanishes is quite contrary 
to our physical expectations inasmuch as a charged parti- 
' cle (of charge — e, position r(i), and velocity v(t) at time 
t), say, orbiting in a plane perpendicular to the magnetic 
field B under its Lorentz force should have an orbital mag- 
netic moment M (= — e/2c[r(t) x v(t)l), where c is the 
speed of light, and with a diamagnetic sign as dictated 
by the Lenz's law (see e.g. [5]). Formally, the vanishing 
of the classical diamagnetic moment follows from the well 
known fact that the canonical partition function involves 
the Hamiltonian for the charged particle (coupled mini- 
mally to the static magnetic field) and a simple shift of the 
canonical momentum variable in the integration makes the 
partition function field-independent, giving zero magnetic 
moment [3]. 
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Physically, the vanishing of the classical diamagnetism 
is due, however, to a subtle role played by the boundary 
of the finite sample [1-4]. It turns out that the diamag- 
netic contribution of the completed cyclotron orbits of the 
charged particles orbiting around the magnetic field in a 
plane perpendicular to it, is cancelled by the paramag- 
netic contribution of the incomplete orbits skipping the 
boundary in the opposite sense in a cuspidal manner. The 
cancellation is exact, and that is the surprise. This cancel- 
lation was demonstrated explicitly some time back [6] for 
the case of a harmonic-potential (V(r) = kr 2 /2) confine- 
ment, which is equivalent to a soft boundary, and finally 
letting the spring constant k go to zero. The treatment 
was based on the classical Langevin equation [7] , and the 
magnetic moment M = — e/2c[r(i) x v(i)] was calculated 
in the infinte-time limit - the Einsteinian approach to sta- 
tistical mechanics. The ordering of the two limits namely 
k — > (the deconfinement limit) and t — > oo (the infinite 
time limit), however, turned out to be crucial and phys- 
ically meaningful - one must let t — > oo first and then 
let k — ► 0. This ensures that the particle is affected by 
the boundary, or the confinement. Thus one had to con- 
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elude that any orbital diamagnctism observed in an ex- 
periment is essentially of quantum-mechanical origin, as 
indeed was derived first by Landau [8]. In the quantum 
case, the above cancellation of the bulk and the bound- 
ary contributions turns out to be incomplete. But again, 
the order of the two limits is all important and was im- 
plicit in the treatment of Landau [3] . It was shown more 
explicitly by Darwin [9]. In fact, one could just use the 
quantum Langevin equation [10] and derive essentially the 
Landau result by properly taking the above 'Darwin limit'. 
The calculated diamagnctic moment is, however, found to 
depend on the frictional term occuring in the quantum 
Langevin equation [11,12]. 

The subtle but essential role of the boundary in all these 
treatments has motivated us to examine the diamagnctism 
for a classical system that has no geometrical boundary - a 
finite unbounded system such as a charged particle moving 
on the surface of a sphere under the appropriate Langevin 
dynamics in the presence of a uniform external magnetic 
field. We were pleasantly surprised to find that the nu- 
merically computed orbital magnetic moment turned out 
to be non-zero, and indeed diamagnctic. To the best of our 
knowledge, this is the first example reported on non-zero 
orbital diamagnctism in principle for a classical system. 
It arises explicitly from avoided cancellation as the sys- 
tem has no boundary. 

Consider a charged particle (charge -e and mass m, an 
electron say) moving on the surface of a sphere of radius 
'a', in the presence of a uniform externally applied mag- 
netic field B directed along the z— axis. The particle mo- 
tion is described by the following classical Langevin equa- 
tion 

</v , . (!) 
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= — (v x B) — Tv + y/2Tk B Tf(t) 
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where L is the friction coefficient, k B T the thermal energy, 
and f is a zero-mean 5— correlated Gaussian random noise, 
i.e., (f a (t)fp(t')) = 5 a p5(t — t'). We recall here that in this 
real space-time (Einsteinian) approach to statistical me- 
chanics, the long-time limit (t — ► oo) of the above stochas- 
tic evolution is expected to give the thermal-equilibrium 
properties. Note that there is no modification of the dis- 
sipation (L) and the related noise term (f(i)) due to the 
magnetic field [7]. 

Specializing now to the spherical-polar coordinates ap- 
propriate to the motion on the surface of the sphere 
(r = a, 0, 0), the Langevin equation reduces to 
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where f , 6 and tf> are the unit vectors directed along the ra- 
dial (r), polar (9) and the azimuthal ((f>) directions. Also, 
f$ and f$ are the forcing noise terms acting along the 9 
and the (f> directions respectively. More conveniently, we 
re-write equation (2) in the dimensionless form 



9 — sin 9 cos ( 



- — sin 9 cos 9(f> — 6 + ^/rjfe , (3) 
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sin 9<p + 2 cos 99(f) = — cos 00 - sin 00 + Jrjfj,, (4) 
7 

where we have introduced the cyclotron frequency lo c = 
eB/mc, the frictional velocity relaxation rate 7 = T/m, 
the thermal forcing strength r\ — 2k B T '/ '{ma 2 "] > ), and the 
dimensionless time r = "ft. Note that r\ is also a dimen- 
sionless quantity. Here overhead dots denote differentia- 
tion with respect to the dimensionless time r. The phys- 
ical quantity of interest is the ensemble averaged orbital 
magnetic moment 

(M(r))=-^ 7 a 2 (sin 2 0(r)^(r)) (5) 

in the long-time limit, where (• • • ) denotes the ensem- 
ble average over the different realizations of the stochastic 
forces fe and fy. 

We now rewrite the above second-order differential 
Langevin equations (4) as four coupled first-order equa- 
tions for 0, x (= 0), (f> 1 and y (= </>), which are then solved 
numerically using a simple Euler-Maruyama scheme [13] 
with a time-step At = 10~ 2 . Averages are evaluated over 
n = 10 6 noise realizations. The number of realizations, 
though quite large, is necessarily finite, and so we resort 
to double average ((•••)) denoting averaging over the en- 
semble as well as over time. This gives for the equilibrium 
magnetic moment 



M. 



eg 



«M(r)» 



1 



r 

Jo 



(M(r)) dr (6) 



as r max -^00. In the context of numerical simulation, 
we have to be careful at the singular polar points = 
and = 7r, where l/sin0 diverges. This is regularized by 

replacing sin0 by v sin 2 9 + e where e is a small positive 
quantity taken to be of order At. Further, inasmuch as 
the physical motion is restricted to < < ir and < 
4> < 2ir, while mathematically however, the equations (4) 
can evolve outside these bounds, we have to set in our 
numerical simulation the following conditions: If 0(t) < 
0, then 0(r) -0(r), x(t) -> -x(t), 0(t) -» <f>(r - 
At) + tt; and if 0(t) > tt, then 0(t) -» 2tt - 0(t), x(t) -» 
—x(t), 4>(t) — > 4>(t — At) — n. This takes care of the 
trajectories that happen to pass through the poles. The 
choice of initial conditions on and (f>, and their time 
derivatives, turns out to be irrelevant for the long-time 
ensemble averaged behavior as indeed is validated by our 
numerical simulation. 

In Fig. 1, we have plotted the dimensionless magnetic 
moment (/x(r)) = 2c/(eja 2 ) (M(t)) as a function of r for 
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certain choice of oj c /j and r\. As can be readily seen, 
the moment is diamagnctic and odd in the magnetic field. 
Also, it can shown to be independent of the sign of the 
charge (electron or hole) as indeed it must be. The fluc- 
tuations seen in the figure are statistical fluctuations due 
to the finiteness of the number of realizations used for en- 
semble averaging. These are thus statistical fluctuations - 
these will, and indeed do, decrease with increasing number 
of noise realizations n. 

For completeness, we have also plotted the computed ve- 
locity distributions (the 9 and cj> components) in fig. 2 and 
these are seen to be essentially Maxwellian as expected, 
with the correct mean square values consistent with the 
fluctuation-dissipation theorem. In fact there is no depen- 
dence on the external magnetic field. 

Fig. 3 shows the variation of the dimensionless magnetic 
moment [i eq (corresponding to M eq ) with the magnetic 
field uj c /-f (which is proportional to B). The plot shows 
an essential linear response which is diamagnetic. 

In Fig. 4, we have plotted the probability density P(/x) 
of the statistical mechanical fluctuations about the equi- 
librium value fi eq . The distribution for the chosen values 
of the parameters is quite broad relative to the mean. (The 
corresponding plot for a system with a boundary is indeed 
known to be broad [14]. Of course, in that case the mean 
is zero.) 

We now return to the main point of this puzzle, namely 
that the classical Langevin dynamics for this finite un- 
bounded system gives a non-zero diamagnetic moment, 
and yet a straightforward calculation using the canonical 
partition function with a minimally coupled Hamiltonian 
gives a free energy which is independent of the magnetic 
field, and therefore, a zero field-derivative of the latter 
implying zero diamagnetism. Now, the classical Langevin 
dynamics provides a real space-time picture of the charged 
particle motion under the influence of fluctuations and the 
concomitant dissipation. Its long-time limit is expected 
to describe thermal equilibrium. And sure enough it does 
reproduce the Maxwellian velocity distribution (fig. 2). 
Moreover, it is manifestly gauge-invariant because it in- 
volves the magnetic field directly without invoking a vector 
potential (indeed, in classical electrodynamics, the vector 
potential is essentially a matter of convenience, unlike in 
the case of quantum mechanics). Also, the computed dia- 
magnetism is consistent with the Lenz's law. On the other 
hand, the canonical treatment based on the Hamiltonian 
underlies all of classical statistical mechanics, but it gives 
zero classical diamagnetism. Our resolution of this puzzle 
is as follows. It is known that classically the static mag- 
netic field does no work on the moving charge inasmuch 
as the Lorentz force (— | (v x B)) acts perpendicular to 
the instantaneous velocity vector. However, such a gy- 
roscopic force can still alter the motion of the charged 
particle so as to give a non-zero magnetic moment without 
changing its energy. In fact, it induces a correlation be- 
tween the velocity and the transverse acceleration due to 
the Lorentz force. Clearly, such a subtle dynamical cor- 



relation, without change of energy, is not captured by the 
equilibrium partition function. But, the Langevin dynam- 
ics manifestly treats this gyroscopic Lorentz force through 
the equation of motion. It is thus our view that the real 
space-time treatment based on Langevin dynamics takes 
into account these subtle correlations involving velocity 
and the transverse acceleration caused by the Lorentz 
force without changing the energy, which is missing from 
the usual partition function. Of course this disagreement 
is only for the special case of a strictly unbounded classical 
system in an external magnetic field. 

Thus, we are forced to admit that there are the fol- 
lowing two alternatives - either there is indeed non-zero 
classical diamagnetism for an unbounded finite system as 
under consideration, or the classical Langevin dynamics 
fails to describe in the long-time limit the thermal equi- 
librium as described by the classical partition function in 
the presence of a magnetic field [15]. Either way, we have 
a non-trivial result. 

Finally, it will be apt at this stage to make a few com- 
ments, bearing on the physical realizability of such a clas- 
sical system. 

First, we recall that the non-zero diamagnetic moment 
for the classical system discussed above is due entirely to 
the absence of a boundary - the avoided cancellation for 
a finite but unbounded system. Now, for the case of the 
quantum mechanical (Landau) diamagnetism too there is 
a cancellation, but it is incomplete [3]. Hence the small- 
ness of the Landau diamagnetism in general. We may 
reasonably expect then that the quantum mechanical dia- 
magnetism for a finite but unbounded system too should be 
different, probably larger, because of the avoided cancella- 
tion [16]. Second, our classical treatment is valid in princi- 
ple for a finite but unbounded system, i.e., a strictly closed 
two-dimensional surface (the charged particle moving on 
the surface of a sphere). It may, however, be physically 
realized to an approximation. Thus, we could consider 
a dielectric microsphere coated with an ultrathin layer of 
a conducting material having small carrier concentration 
at room temperature, e.g., a non-degenerate system with 
the degeneracy temperature much smaller than the room 
temperature, as in the case of a doped high-mobility semi- 
conductor. (By ultrathin we mean here a thickness <C the 
thermal dc Broglie wavelength of charge carriers so as to 
freeze out the radial motion quantum mechanically, mak- 
ing the system essentially a two-dimensional classical gas 
of charged particles moving on the surface of the sphere. 
Thus, quantum mechanics helps us realize a dimensional 
reduction - an essentially two-dimensional closed surface). 
We could then consider a finite volume fraction of an in- 
ert medium (paraffin say) occupied by the above micro- 
spheres. This system should have a measurable diamag- 
netic response which will be essentially classical. Third, 
it should be interesting to consider more general geome- 
tries such as that of a triaxial ellipsoid where different axes 
ratios can mimick very different physical situations. Per- 
haps much more interesting will be to try out topologies 
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Fig. 1: (Color online) Plot of the ensemble averaged dimension- 
less magnetic moment (j.(t) as a function of the dimensionlesss 
time t for uJc/'y = ±10.0 and n = 1.0. Clearly the moment can 
be seen to be odd in the magnetic field B and is diamagnetic. 
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Fig. 2: (Color online) The velocity distribution on a sphere 
in the presence of magnetic field computed from the Langevin 
dynamics simulation. The dashed curve is the corresponding 
Maxwellian distribution. Here io c j^ — 10.0 and 77 = 1.0. 



other than that of a sphere and look for qualitative dif- 
ferences [16]. The numerically estimated value of the the 
diamagnetic moment for a charged particle (say, electron) 
moving on a sphere of radius a = 100/im for B ~ 5 kilo- 
gauss, 7 ~ 10 9 s _1 turns out to be ~ 1 Bohr magneton 
per electron which is quite large. Thus we may have a gi- 
ant classical diamagnetism. Of course, the measured bulk 
susceptibility for the physical classical system suggested 
above may have much smaller values because of the re- 
alizable parameter values. But, the point of principle at 
issue will have been made. It is our hope that experimen- 
talists may take note of this possibility. 

To conclude, we have shown that a classical system of 
charged particles moving on a finite but unbounded surface 
(of a sphere), as described by Langevin dynamics, has a 
non-zero orbital diamagnetic moment. This moment can 
be large. 
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Fig. 3: (Color online) Plot of the dimensionless magnetic mo- 
ment fi eq as a function of uj c /"/ (proportional to the magnetic 
field) for 7] = 1.0. Again the moment is found to be odd in B 
and is diamagnetic in sign. 
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